Simulation of Grover's quantum search algorithm in 
a Ising nuclear spin chain quantum computer with 
first and second nearest neighbour couplings 

G V Lopez, T Gorin and L Lara 

Departamento de Fisica, Universidad de Guadalajara 
Blvd. Marcclino Garcia Barragan y Calzada Olimpica 
4480 Guadalajara. Jalisco, Mexico 

E-mail: gorin@pks.nipg.de 

Abstract. We implement Grover's quantum search algorithm on a nuclear spin chain 
quantum computer, taking into Ising type interactions between nearest and second 
nearest neighbours into account. The performance of the realisation of the algorithm 
is studied by numerical simulations with four spins. We determine the temporal 
behaviour of the fidelity during the algorithm, and we compute the final fidelity as 
a function of the Rabi frequency. For the latter, we obtained pronounced maxima at 
frequencies which fulfil the condition of the 27rfc-method with respect to the second 
nearest neighbour interactions. 



PACS numbers: 03.67.-a, 03.67.Lx 



Simulation of Grover's quantum search algorithm 



2 



1. Introduction 

The Ising-spin chain quantum computer is a simple theoretical model system, which has 
many of the relevant features of a physical realisation of a quantum computer. The bare 
model is defined in terms of a time-independent Hamiltonian Hq, describing a system 
of spins in a strong magnetic field and coupled by the Ising-interaction. Additional 
time-dependent perturbations {e.g. radio- frequency pulses, which may be switched on 
and off at will) are used to implement certain elementary quantum gates, from which 
any unitary quantum algorithm can be constructed (universal quantum gates). 

Originally, this model has been proposed in P, El [3l H], and a large number of 
theoretical studies has been devoted to it (to mention a few important works: [5l[6l[7]). 
As a physical realisation, one may think of nuclear spins embedded in a solid state system 
in an external strongly inhomogeneous magnetic field. Ultracold atoms in optical lattices 
may provide a possibly less demanding alternative [8] . 

In the present work, we assume an Ising-interaction between nearest and next- 
nearest (second) neighbours [9l[T0]. In this case, the bare Hamiltonian becomes 

n n— 1 2 

-Ho = -j2'^k n -2jj2 n - 2 / n 1^2 , (i) 

fc=l k=l k=l 

where Wk is the Larmor frequency of spin k. Since the dipole interactions may be 
induced by electron clouds around the nuclear spins, it is realistic to consider J and J' as 
independent parameters. The Hamiltonian Hq is diagonal in the so called computational 
basis, which consists of products of single spin (qubit) eigenstates of J|, which are 
denoted by |Ofc) and |lfc) respectively: 

II Wk) = \a,) . (2) 

We denote the 2" eigenstates of Hq by ■ ■ ■ , ao)} with ak = 0,1 for k = 

0,---,n- 1. 

In order to implement the desired quantum gates, we use monochromatic 
electromagnetic radio-frequency pulses (RF-pulses) described by the term W{t), 

H = Ho + W{t) , W{t) = \Y. (e'^"*+^^4+ + e-^('^*+^)4-) , (3) 

k=o 

where the frequency the phase offset and the Rabi frequency VL are free parameters 
(n/7 is the amplitude of the electromagnetic field and 7 is the gyromagnetic ratio for the 
spins). = {IkV = |Ofc)(lfc| is the raising operator of qubit k. The desired quantum 
gates (and the whole quantum algorithm) are realised within the interaction picture. We 
apply the above defined RF-pulses with a rectangular envelope, choosing the frequency 
w in resonance with an appropriate transition in Hq. The resulting unitary evolution is 
denoted by 

i?r((^, ^) , w = Wk + fiJ + iyJ' , (4) 
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where /i and u are integers. The angle (p gives the phase offset of the RF-pulse, whereas 
'd = Qt determines the duration r of the pulse. The RF-pulses are selected in such 
a way that they yield the desired quantum gates exactly, if all off-resonant transitions 
are neglected (resonant approximation). For the "real" evolution will be obtained from 
numerical simulations, where all transitions are taken into account. The deviation from 
the "ideal" evolution is quantified by fidelity [11], the absolute value squared of the 
wave-function overlap for both evolutions. In order to minimise the decay of fidelity. 
We apply the 27r/c-method to suppress those transitions which are nearest to resonance. 
These are of order J', and therefore related to the second neighbour interaction [TO] . 

In the present work, we study the effects of unitary errors (due to non-resonant 
transitions) on a quantum algorithm of intermediate length. Such an algorithm may 
consist of several hundreds of pulses, even though it is still far away from the regime, 
where a quantum computer could really demonstrate its superiority over a classical one 
(after all, our work relies on numerical simulations on a classical PC). It is however also 
far beyond pulse sequences, which only realise one or a few elementary qubit gates. Of 
particular interest will be the question whether we can from the fidelity decay of few 
gate operations extrapolate to the fidelity decay of the whole quantum algorithm. 

As an example we consider Grover's search algorithm [121 [13] which has also been 
implemented experimentally [T^ [T5l [T6] . We will implement this algorithm on a spin 
chain of length four, where one of the spins will serve as the "ancilla" qubit. The 
basic elements of this algorithm are single-qubit Hadamard gates. Due to the Ising 
interactions, these gates are quite expensive in the Ising-spin chain quantum computer 
as we will see below. The big majority of the pulses are used to implement these gates. 
The other gates are n-qubit phase-gates, where n is the number of bits of the search 
space (in our case n = 3). 

In the following section (section [2]) , we will discuss Grover's search algorithm as 
such, and its implementation in the Ising-spin chain quantum computer. In section [3l 
we show numerical simulations of the search algorithm, discuss the fidelity decay during 
the execution of the algorithm, and study this fidelity as a function of Rabi frequency. 
Conclusions are provided in section HI 

2. Realisation of Grover's search algorithm 

In its most basic form [T7], Grover's quantum search algorithm requires n qubits 
to prepare the "enquiry" states (in the data register) subsequently presented to the 
"oracle", and a single additional qubit (the "ancilla" qubit) required by the oracle to 
communicate its answer. Starting from the ground state of the system, |0„_i, ■ ■ ■ , Oq), 
the algorithm itself consists of an initial step, where the complete superposition of all 
basis states is created by applying Hadamard gates to all qubits of the data register. 
Let us denote a single Hadamard gate applied to qubit j by Hj and the sequence of 
Hadamard gates to create the superposition state in m qubits by H®^. The creation of 
the superposition state is followed by a Grover operator, which may be repeated several 
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times, until the data register contains the desired target state with sufficiently high 
probability. 

The Grover operator is again composed of Hadamard gates and conditional phase 
reflections: 

G = 0{a) H^"" So if®" , (5) 

where a may be any number representable in the data register. The unitary oracle 
0(a), leaves all basis states untouched, except for the state \a), to which it applies a 
sign change (phase reflection). The other unitary operator 5*0, in turn, changes the sign 
of all basis states but leaves the state |0) untouched. While the Hadamard gates are 
single qubit gates, 0{a) and 5*0 are true ra-qubit gates. 

In this work we implement the Grover algorithm in a quantum register of four 
qubits: {^3, "^i, "^o)}- For numerical simulations, it is convenient to use the qubits 
/c = 0, 2, 3 as data register, and the qubit k = 1 as the ancilla qubit. In this way, the 
ancilla qubit is coupled directly (by first or second neighbour couplings) to all qubits 
in the data register. In what follows, we discuss the implementation of the Hadamard 
gates, the conditional phase reflection 5*0 and the oracle, in separate subsections. 



2.1. Hadamard gates 

The Hadamard gate is a single qubit gate. In the Hilbert space of that qubit the gate 
acts unitary matrix: 

1 1 
1 -1 




(6) 



It can be decomposed into elementary qubit rotations as follows: 

H = Ry{7r) Ry{TT/2) R,{7r) , (7) 

where R^{'d) = R{7t,^), Ry{^) = R{'k/2,^) and 

cos(^/2) ie^^ sin(^/2) \ 
^^^'^^-|^ie-'-sin(^/2) cos(^/2) j' 

The single qubit rotations of the form (IHl) may be obtained by applying appropriate 
electromagnetic pulses as defined in ([3]) and (jl]) to the Ising spin chain. In the presence of 
neighbouring qubits, one has to repeat that pulse sequence for all possible configurations 
of the neighbouring qubits, in order to obtain the true single qubit gate, independent 
of the state of the neighbouring qubits. Denoting the Hadamard gate applied to qubit 
i as Hi., we obtain specifically: 

Ho= w n Ki^/'^.^) Ki^/'^.^m (9) 

lJ,= -l,l u = -l,l 

H2= W n ^r(vr/2,7r)i?f(7r/2,7r/2)i?f(7r,7r) (10) 

Ai=-2,0,2 u=-l,l 

^3= n n i?r(^/2,vr)i?r(^/2,vr/2)i?r(vr,vr), (11) 
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which implements the unitary gate if®'^ = HqH2H-^ on the data register. 

2.2. The conditional phase inversion Sq 
Formally, the gate 5*0 may be described by 

5o|a) = (-l + 25,,o)|«). (12) 

The (unconditional) phase inversion of the state of a single qubit may be obtained by 
applying the pulse -R(0, 27i) to it. In the present case, we obtain 5*0 by applying i?(0, 2%) 
to the qubit i = 1 repeatedly - covering all possible configurations of its neighbouring 
qubits except for that case, where the data-qubits form the state |0). This leads to the 
pulse sequence 

So = rI'-\(},2it) i??'^(0,27r) i??'~^(0, 27r) i?^^'~^(0, 27r) Rf'\0,2iT) . (13) 

2.3. The oracle 

The oracle is a unitary gate, which changes the sign of the state of the data register, 
depending on its contents. In the search algorithm, it is thought that any number |a) 
representable in the data register serves as an index to given items in a data base. The 
search consists in finding the index of that item which has a certain unique property 
{e.g. the telephone number of a certain person with a certain address). In the present 
case, we will request that the index itself be a certain number out of the set {0, 5, 8, 13} 
(decimal notation). In those cases, the oracle may be implemented by a single 27r-pulses: 



O{0) = RI'\0,27t) (14) 

0(8) = i?i'"^(0,27r) (15) 

0{5) = Rf'\0,27r) (16) 

0(13) = i?i~^'"^(0,27r) . (17) 



The other numbers representable in the data register are {1,4,9,12}, but a oracle 
selecting those states would require many more qubits. 

3. Numerical simulations 

Our simulations normally start with the ground state, |0000). Then, we generate the 
total superposition state in the data register and perform two Grover steps. This is 
sufficient to obtain the desired target state with very high probability (above 95%). With 
the pulse sequences given in section [2], one needs 42 pulses for the initial superposition 
state and 90 pulses for each Grover operator step. 

Unfortunately, the pulse sequence described in section [2] is not optimal in terms of 
accuracy. The direct implementation of the Hadamard gates according to the scheme 
described there, leads to prohibitively large errors. Luckily, these errors can be reduced 
efficiently, by recombining the pulses in the sequence for one Hadamard gate into 
shorter sequences implementing several single qubit gates. This can be done due to 
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Figure 1. The absolute value squared of the fidelity during the Grover algorithm with 
two Grover operator steps. Time is measured in units of the duration of a single 7r/2- 
pulse. The algorithm is simulated for different target states \a) (see legend). The Rabi 
frequency is f2 = 0.1008, which suppresses non-resonant transitions as far as possible. 



the commutativity of the operators appearing in Thus, for e.g. the Hadamard gate 
Hq, we obtain the following pulse sequence 

Ho = {U,,.=-,,iKi^/'2,n)] {n,,=_i,ii?rW2,vr/2)} {n,,.=-i,i (^r, vr)} , (18) 

with much smaller errors. We use similar rearrangements for H2 and H^. 

The overall error can be further decreased by rearranging the outer sequences of 
TT-pulses in a similar manner into sequences of single qubit 7r/2-pulses. These measures 
drastically increase the number of pulses for the implementation of the algorithm. Now, 
we need 20 pulses for the Hq and the Hs gate and 30 pulses for the H2 gate. In total 
this gives 70 pulses for the H'^^ gate and 146 pulses for each Grover operator step. Note 
however that the amount of time to implement those gates is the same as before. 

In this section, we are mainly concerned with unitary errors due to non-resonant 
transitions, and due to phase-errors which are not suppressed by the 27r/c-method. 
Thus, we will study in some detail the behaviour of fidelity during the algorithm. The 
fidelity [11], f(t), is the overlap between the state evolving under the ideal evolution 
(where only resonant transitions are taken into account) and the state evolving under the 
true evolution governed by the Hamiltonian in ([3]). The fidelity (or fidelity amplitude) 
is a complex quantity. Thus, we will normally study its absolute value squared, 

Figure [H shows during the whole search algorithm for four different target 

states, selected by different oracle pulses: red curve (a = 0), blue curve (a = 5), light 
blue curve {a = 8), and pink curve {a = 13). In these simulations, we set the Rabi 
frequency to ^2 = 0.1008, which corresponds to the 27r/c-method with k = 4. Since k is 
even, non-resonant transitions are suppressed during vr-pulses, as well as 7r/2-pulses. 
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Figure 2. The absolute value squared of the fidelity during the Grover algorithm 
with two Grover operator steps. Time is measured in units of the duration of a single 
7r/2-pulsc. The algorithm is simulated for different target states \a) (see legend). The 
Rabi frequency is = 0.135, which yields similarly high fidelities as in figure [TJ The 
smooth blue line shows an exponential, fitted to the initial decay (V'''7r/2 < 70). The 
fastest decay, roughly following the exponential is obtained for random initial states 
(lowest red curve). 

Figure [2] shows for the same Grover algorithm, but for the Rabi frequency, 

Q = 0.135, which corresponds to the 27r/i;-method with k = 3. Since k is odd in this 
case, one would expect that non-resonant transitions are suppressed during vr-pulses, 
but not during 7r/2-pulses. Still, we find similarly high fidelities as in the previous case. 
This surprising result is discussed in more detail, below. In this figure, the fidelity along 
the full algorithm does not follow the initial trend of relatively fast decay. To obtain 
that initial trend, we fitted an exponential to the initial preparation step (the first 70 
7r/2-pulses). 

This non-generic behaviour is interesting, because it shows that it is not sufficient to 
characterise the performance of individual quantum gates by their average fidelity-loss 
rates. With no other information than such loss rates, one would greatly underestimate 
the overall fidelity of the algorithm. A convenient way to obtain average loss rates, is to 
use random initial states. The lower red curve shows averaged over 100 different 

initial states chosen at random from the invariant (under unitary transformations) 
ensemble of normalised states. The curve still shows large fluctuations which should 
disappear in the limit of a larger sample size. However, the general trend of a global 
exponential decay is clearly observable. It also agrees quite nicely with the exponential 
flt (smooth blue curve). 

In flgure [3] the Rabi frequency is chosen to allow non-resonant transitions, Q = 
0.2550. By consequence, we flnd a much faster decay of the fldelity. In distinction to 
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Figure 3. The fidelity (its absolute value squared) during tlie Grover algorithm with 
two Grover operator steps. Time is measured in units of the duration of a 7r/2-pulse. 
The algorithm is simulated for different target states \a) (see legend). 
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Figure 4. The fidelity (blue line) as well as its absolute value squared (red line) of 
the whole Grover algorithm as a function of the Rabi frequency, Q. The algorithm is 
simulated for the target state |a), a = 13. 

the previous figures, the fidelity now follows much closer an exponential decay, although 
there are still very large fiuctuations between the different cases (different oracles). 

In figure. H] we plot the real part of the fidelity (blue line) as well as its absolute 
value squared (red line) for the whole algorithm with two Grover operator steps as a 
function of the Rabi frequency, Q. We chose the oracle 0(13), which selects the target 
state I a) with a = 13. One can clearly see pronounced peeks of maximal fidelity, where 
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the 27r/c-method is at work. The Rabi frequencies obtained from the condition of the 
27r/c- method are indicated by vertical hnes, where k increases from k = 1 near Q ^ 0.45 
up to /c = 6 near Q ^ 0.067. 

For integer k, the 27r/c-method suppresses non-resonant transitions only for yr-pulses. 
In order to suppress non- resonant transitions also for 7r/2-pulses, k must be even. Since 
the present algorithm consists mainly of 7r/2-pulses, we should expect rather small 
fidelities at Rabi frequencies which correspond to odd values for k. Nevertheless, even 
for odd values of k, we find pronounced maxima of maximal fidelity. They are almost as 
high as the maxima corresponding to even fc's. Somehow, it happens that non-resonant 
transitions which occur during one 7r/2-pulse are undone during a second 7r/2-pulse. 
While this might not be so surprising in the case of the vr-pulses which have been split 
deliberately into 7r/2-pulses, it is more surprising in the case of Ry{7r/2) which also forms 
a part of the Hadamard gate [see equation (^^], and which is implemented with "true" 
7r/2-pulses. 

4. Conclusion 

We implemented Grover's quantum search algorithm on a Ising spin chain of length 
four, taking into account first and second neighbour couplings. The coupling strength 
of second neighbours is much weaker than that of first neighbours, which means that 
the dominant source for errors are near resonant transitions, where the detuning is of 
the order of J'. In other words, the second neighbour couplings are the dominant source 
for errors in our setup. The accumulated error in the course of the quantum algorithm 
are quantified by the fidelity. We investigated its temporal behaviour as well as its final 
value, at the end of the algorithm, as a function of the Rabi frequency. 

We found that by properly rearranging the RF-pulses, the errors in the 
implementation of the Hadamard gates are greatly reduced. By properly choosing the 
Rabi frequency we obtained fidelities beyond 90%, for the whole algorithm with two 
Grover steps. This corresponds to an average loss of fidelity of the order of 4 x 10~^ (per 
TT-pulse) and 6 x 10~^ (per Hadamard gate), numbers which are compatible with recent 
requirements for fault tolerant quantum computation [18]. However, we also found 
that the loss of fidelity in different subsections of the algorithm varies considerably 
(sometimes even in sign), which means that it is in general not possible to estimate the 
fidelity of the whole algorithm, knowing the fidelities of its parts. 

In the present work, we implemented the 27rA;-method to suppress non-resonant 
transitions. Since the big majority of the RF-pulses used are 7r/2-pulses, we would have 
expected high fidelities at Rabi frequencies corresponding to even k, but rather low 
fidelities at Rabi frequencies corresponding to odd k. Surprisingly, we found similarly 
pronounced maxima of the fidelity at Rabi frequencies corresponding to even and odd 
values for k. 

Preliminary studies using fluctuating Larmor frequencies (simulating tiny variations 
of the external magnetic field) indicate that these findings are not sensitive to external 
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noise. A more detailed study of the effect of noise will be the subject of future 
investigations. 
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